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HIGHLIGHTS 


•  40  Ah,  Saft  batteries  were  studied  using  EIS  at  0%  SOC  as  function  of  temperature. 

•  The  batteries  were  modelled  using  the  equivalent  circuit  approach. 

•  EIS  spectra  exhibit  superimposed  impedance  arcs  that  are  temperature  dependent. 

•  Non-viable  and  viable  batteries  could  be  discriminated. 

•  Activation  energy,  polarisation  resistance  and  time  constant  could  be  extracted. 
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A  theoretical  equivalent  circuit  model  of  two  40  Ah  commercial  batteries  supplied  by  Saft  was  formu¬ 
lated  to  interpret  electrochemical  impedance  spectra  as  a  function  of  temperature  at  zero  state-of- 
charge.  The  batteries  were  chosen  to  represent  viable  and  non-viable  products  from  the  Saft  produc¬ 
tion  line.  This  paper  is  the  first  contribution  from  a  project  in  Newcastle  to  develop  a  rapid,  non¬ 
destructive  analytical  method  for  quality  control  on  the  battery  production  line.  Using  the  model,  it 
proved  straightforward  to  discriminate  between  viable  and  non-viable  batteries  based  on  the  temper¬ 
ature  dependence  of  the  electrochemical  reaction  kinetics  at  the  electrodes,  the  average  diffusion  co¬ 
efficients  and  the  charge  transfer  resistances.  Furthermore,  as  well  as  marked  differences,  for  example, 
between  the  time  constants,  activation  energies,  polarisation  resistances  of  the  viable  and  non-viable 
batteries,  the  individual  contributions  of  the  anodes  and  cathodes  in  the  batteries  were  de-convoluted 
and  interpreted  through  the  model  framework. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Secondary  lithium  ion  batteries  are  ubiquitous  in  portable 
electronic  devices  and  are  finding  increasing  application  in 
hybrid  electric  vehicles  [1,2].  During  the  charging  of  such  devices, 
Li  ions  de-insert  from  the  cathode  [3]  (typically  comprised  of 
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metal  oxides  of  the  form  LixM02)  and  are  inserted  into  graphite 
anodes  [4]  and  vice  versa  on  discharging.  Due  to  the  very  high 
energy  density  of  lithium  ion  batteries,  quality  control  is  of 
paramount  importance  in  the  production  of  these  devices  to 
ensure  safety  as  well  as  durability.  In  order  for  these  re¬ 
quirements  to  be  met,  it  is  essential  that  faults  in  batteries  due  to 
manufacturing  errors,  etc.  can  be  detected  rapidly  and  with  high 
accuracy  (i.e.  with  low  incorrect  rejection  rates);  preferably  on 
the  production  line. 

Batteries  are  generally  characterized  by  using  standard  cur¬ 
rent/voltage  (I IV)  measurements  [5]  and  by  the  more  sophisti¬ 
cated  Electrochemical  Impedance  Spectroscopy  (EIS)  [6],  The 
latter  involves  fitting  EIS  data  to  a  theoretical  model,  based  on 
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Fig.  1.  Schematic  representation  of  electrified  interface. 


circuitry  comprising  electronic  components  which  may  be 
correlated  with  specific  elements  (e.g.  anode,  cathode,  electro¬ 
lyte,  etc.)  and  processes  (e.g.  diffusion)  in  the  battery  [7,8].  In 
this  way,  resistance  due  to  the  Solid  Electrolyte  Interface 
(SEI)  layer,  the  charge  transfer  processes  of  the  anode  and 
cathode,  the  electrolyte  and  connector  impedances,  diffusion 
effects  [9,10]  and  changes  due  to  the  electrochemical  and 
chemical  reactions  [11]  in  the  battery,  etc.  (see  Fig.  1)  can  be 
deduced  from  a  comparison  of  the  theoretical  model  with 
experimental  data.  The  rate  constant  of  fast  electron  reactions 
can  be  also  monitored  [12], 

EIS  is  also  called  AC  impedance  spectroscopy  [13]  and  was  first 
applied  in  the  field  of  wet  electrochemistry.  It  has  been  employed 
in  corrosion  studies  where  complicated  processes  involve  surface 
and  solution  reactions.  Moreover  the  technique  has  also  been  used 
extensively  for  lead  acid  battery  characterisation.  EIS  often 


Table  1 

List  of  symbols. 

Symboi/unit 

i/H 

Inductance 

R/Q 

Resistance 

Rp/Q 

Polarisation  resistance 

Ret /Q 

Charge  transfer  resistance 

R„/U 

Warburg  resistance 

CPE  (Q) 

Constant  phase  element  (capacitance) 

r/s 

Time  constant 

Tw/S 

Warburg  time  constant 

T/K  or  Celsius 

Temperature 

Vo/S  s1'2 

Admittance 

w/ rad  s-1 

Angular  frequency 

I< 

Rate  constant 

A 

Pre-exponential  factor 

Ea/J  mol  1 

Activation  energy 

R/J  moD1  IC1 

Universal  gas  constant 

cr/Q  s-1 

Diffusion  coefficient 

zw 

Warburg  impedance 

SOC 

State  of  Charge 

Table  2 

Electrical  characteristics  of  the  supplied  batteries. 


Nominal  voltage 

3.6  V 

Typical  capacity  (C  rate,  +20  °C) 

41  Ah 

Minimum  capacity  (C  rate,  +20  °C) 

40  Ah 

Maximum  voltage  (continuous) 

4.0  V 

Maximum  voltage  (pulse) 

4.1  V 

Minimum  voltage  (continuous) 

2.7  V 

Minimum  voltage  (pulse) 

2.0  V 

becomes  a  technique  of  choice  to  study  the  electrified  interface 
and  hence  has  been  used  in  lithium  ion  batteries.  Briefly,  the  EIS 
literature  on  lithium  ion  batteries  may  be  represented  by  typical 
papers  such  as  those  by  Zhuang  et  al.  [14  and  Chagnes  and 
Swiatowska  [15],  Chang  and  Park  [16]  have  described  and  dis¬ 
cussed  the  development  of  new  EIS  methods  to  reduce  measure¬ 
ment  time. 

Typical  lithium  ion  battery  modelling  using  EIS  is  focused  on  the 
electrolyte  and  the  porous  electrodes/electrolyte  interfaces  in  order 
to  determine  the  optimum  parameters  for  better  batteries  design 
and  to  characterise  ageing  effect  [17],  Valuable  information  for 
diagnosing  faults  whilst  in  operation  may  be  obtained  and  accurate 
State  of  Charge  (SOC)  determined  [18].  There  are  many  models 
utilising  equivalent  circuits  to  describe  the  electrochemical 
behaviour  of  lithium  ion  batteries,  the  simplest  is  the  single  parti¬ 
cles  (SP)  model  [19].  Other  models  include,  in  order  of  increasing 
sophistication  (i.e.  decreasing  maximum  error  rate:  the  Internal 
Resistance  (IR)  model,  One  Time  Constant  (OTC),  Two  Time  Con¬ 
stants  (TTC))  to  try  to  better  monitor  the  cell  dynamics  [20,21  ],  the 
Rint,  the  RC,  the  Thevenin  [22],  the  partnership  for  a  new  genera¬ 
tion  of  vehicle  (PNGV)  [23]  and  the  improved  Thevenin  or  Dual 
Polarisation  (DP)  model  [24],  Of  all  the  models  above,  DP  model 
characterises  most  accurately  the  dynamics  of  the  lithium  ion 
battery  and  provides  the  most  accurate  estimation  of  SOC.  It  was 
found  that  we  could  only  fit  the  data  from  the  Saft  batteries  with 
multiple  time  constants  and  hence  we  employed  a  variation  of  the 
DP  model. 

Another  consideration  in  the  field  of  failure  diagnosis  is  to 
monitor  the  temperature.  Temperature  imbalance  has  been 
observed  in  batteries,  and  this  can  lead  to  unexpected  performance 
and  life  degradation  [25].  Such  problems  are  related  to  battery 
design.  Thus,  in  a  cylindrical  battery  [26],  the  centre  is  slightly 
warmer  than  the  exterior  which  leads  to  preferential  reactions  at 
the  centre  upon  charge  and  discharge.  In  a  prismatic  cell  [27]  the 
temperature  gradient  is  different  due  to  the  difference  in  shape  and 
in  particular  the  fact  that  the  latter  are  generally  very  thin  and 
dissipate  heat  more  effectively. 

In  general,  physical  processes  taking  place  in  galvanic  cells  may 
be  distinguished  on  the  basis  of  their  time  constants,  and  hence 
frequency  dependent  responses.  Thus,  high  frequency  behaviour  is 
dominated  by  the  inductive  response  of  the  connecting  leads,  the 
medium  frequency  range  by  the  charge  transfer  at  the  electrodes/ 
electrolyte  interfaces  and  the  low  frequency  region  by  diffusion  of 
lithium  ions  through  the  electrode  matrix  and  SEI  layer  [28], 

This  paper  presents  initial  studies  aimed  at  developing  a  non- 
invasive,  ElS-based  system  for  the  rapid  characterization  of  vehic¬ 
ular  lithium  ion  batteries.  The  work  reported  in  this  paper  was 


Table  3 

Mechanical  characteristics  of  the  VL41 M  batteries. 


VL41 M  sleeved 

VL41 M  formed 

Diameter 

55.0  mm  max. 

54.6  mm  max. 

Height 

223  mm  max. 

222.8  mm  max. 

Weight 

1.1  kg  max. 

1.1  kg  max. 
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Fig.  2.  Nyquist  plots  of  battery  1  from  EIS  at  0%  SOC  at  15  °C  (  •■■  (>••),  20  °C  (  -A- -■), 
25  °C  (■■■□■■■)  and  30  °C  (•••©••■). 

carried  out  at  zero  state-of-charge  (0%  SOC)  i.e.  3.2  V,  where  SOC  is 
defined  as  the  ratio  of  the  energy  available  to  the  total  energy  ca¬ 
pacity  of  the  battery  [29]  expressed  as  a  percentage,  as  a  function  of 
temperature.  A  subsequent  paper  will  report  studies  as  a  function 
of  SOC  at  constant  temperature.  Two  batteries  (see  Table  1 )  were 
supplied  by  Saft  (Bordeaux  France),  one  chosen  as  representative  of 
a  viable  battery  (battery  2),  the  other  (battery  1 )  chosen  on  the  basis 
of  having  a  significantly  faster  self-discharge  rate  than  the  other 
(i.e.  3.940  mV  day-1  in  comparison  to  2.108  mV  day-1).  Given  that 
both  cells  exhibited  self-discharge,  0%  SOC  was  chosen  as  the  ca¬ 
pacity  at  which  measurements  would  be  taken  in  order  to  provide  a 
fixed  point  of  reference. 

Current  practice  is  for  batteries  to  be  rejected  on  the  basis  of 
enhanced  self-discharge  rates.  This  requires,  typically,  two  weeks  of 
monitoring  compared  to  less  than  10  min  with  the  present  tech¬ 
nique.  With  further  optimisation  and  refining  we  believe  a  test  time 
of  ca.  1  min  is  achievable. 

2.  Experimental  methods 

A  set  of  impedance  spectra  were  recorded  on  the  two  cylindri¬ 
cal,  Saft  Li-ion  VL41M  batteries  at  0%  SOC  and  at  various  temper¬ 
atures.  The  electrical  specifications  of  the  batteries  are  summarised 
in  Table  2  and  the  mechanical  properties  in  Table  3.  Battery  testing 
was  performed  using  a  Bio-logic  HP1005/VMP3B-100  electro¬ 
chemical  impedance  analyser. 

The  batteries  were  placed  in  an  environmental  simulation 
chamber  (Binder  GmbH)  capable  of  maintaining  a  constant 


Fig.  3.  Nyquist  plots  of  battery  2  from  EIS  at  0%  SOC  at  15  °C  (••■O--).  20  °C  (•••A-  -). 
25  °C  and  30  °C  (-O-). 
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Fig.  4.  A  schematic  of  the  equivalent  circuit  for  non-linear  least  square  fitting. 

temperature  over  the  range  -40  to  180  °C.  The  batteries  were 
allowed  to  reach  thermal  equilibrium  over  12  h  before  testing.  The 
primary  aim  of  the  work  reported  in  this  paper  was  to  determine  if, 
and  to  what  extent,  the  EIS  responses  of  the  two  batteries  could  be 
discriminated.  Thermocouples  were  connected  to  the  top  and  the 
bottom  of  each  battery  to  monitor  temperature  variation  within  the 
batteries.  EIS  sweeps  were  performed  at  stable  battery  tempera¬ 
tures  of  15,  20,  25  and  30  °C. 

The  batteries  were  discharged  using  constant  current  (CC)  at  a 
rate  of  C/10  to  the  desired  voltage  3.2  V  followed  by  constant 
voltage  (CV)  until  the  current  fell  to  zero,  to  ensure  equilibrium  at 
0%  SOC.  The  batteries  were  held  at  this  voltage  for  3  h  before  each 
EIS  experiment  was  performed. 

The  EIS  experimental  conditions  were:  perturbation  voltage 
amplitude,  Va  =  3  mV  over  the  frequency  range  1  kHz  to  100  mHz, 
at  ten  points  per  decade. 


Table  4 

Fitting  parameter  values  for  impedance  spectra  measured  at  15,  20,  25  and  30  °C. 


T=  15  °C 

T=20°C 

T=  25  °C 

T=30°C 

(a)  Parameters  for  battery  1 

1,/H  (xl0~6) 

0.2382 

0.1570 

0.3130 

0.2266 

(±0.0061) 

(±0.0077) 

(±0.0059) 

(±0.0087) 

r,/q  (xio-3) 

1.407 

1.397 

1.128 

1.291 

(±0.007) 

(±0.007) 

(±0.004) 

(±0.005) 

CPEa/F  s-"(xl0-3) 

20.91 

8.78 

77.34 

96.52 

(±0.23) 

(±0.10) 

(±0.11) 

(±0.01) 

CPEc/F  s-"(xl0-3) 

1.70 

1.721 

64.27 

7.46 

(±0.27) 

(±0.17) 

(±0.01) 

(±0.04) 

Rp.a/Q(xl0-6) 

1815 

1541 

56.73 

0.386 

(±7.868) 

(±7.007) 

(±3.565) 

(±0.087) 

Rp.c/Q(xl0-6) 

240 

132 

140 

4.031 

(±2.542) 

(±1.813) 

(±1.635) 

(±0.291) 

szw(xio-3) 

5.19 

5.70 

6.61 

7.35 

(±0.05) 

(±0.04) 

(±0.03) 

(±0.04) 

n2  =  n3 

0.86 

0.86 

0.86 

0.86 

X2(xl0-3) 

3.145 

0.903 

0.859 

0.523 

Residual  standard  4.916 

error  (xlCT6) 

(b)  Parameters  for  battery  2 

4.807 

3.27 

3.132 

I,/H  (xl0~6) 

0.2429 

0.2280 

0.2872 

0.2378 

(±0.0058) 

(±0.0061) 

(±0.0171) 

(±0.0054) 

r,/q  (xio~3) 

1.375 

1.325 

1.245 

1.229 

(±0.005) 

(±0.005) 

(±0.008) 

(±0.003) 

CPEa/F  s~"  (xl0~2) 

2.00 

1.85 

1.85 

1.98 

(±0.43) 

(±0.18) 

(±0.04) 

(±0.05) 

CPEc/F  s-"(xl0-2) 

1.43 

1.63 

1.61 

8.53 

(±0.01) 

(±0.01) 

(±0.01) 

(±0.01) 

Rp.a/£1  (xlO-6) 

2543 

1588 

1012 

2.653 

(±7.868) 

(±4.064) 

(±5.657) 

(±2.403) 

Rp.c/Q(xlO-6) 

258 

199 

172 

22.46 

(±7.868) 

(±7.401) 

(±4.664) 

(±0.507) 

szw(xio-3) 

4.80 

5.38 

6.19 

6.99 

(±0.05) 

(±0.08) 

(±0.03) 

(±0.04) 

n2  =  n3 

0.86 

0.86 

0.86 

0.86 

X2(xl0-3) 

1.027 

1.729 

0.647 

0.404 

Residual  standard 
error  (x  10  6) 

5.901 

4.403 

3.176 

3.1 

The  figures  in  brackets  are  the  relative  standard  errors  in  each  parameter. 
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Fig.  5.  Arrhenius  plot  battery  1  ( ■  a  )  and  battery  2  (•--■)  at  0%  SOC. 


The  fitting  of  the  data  so  obtained  was  conducted  using  EC-Lab 
software  (Zfit).  The  computer  program  employed  for  analysis  by  Zfit 
was  based  on  a  complex,  Non-Linear  Least  Squares  (NLLS)  fitting  of 
the  model  parameters  of  the  impedance  spectrum  to  the  measured 
spectrum.  NLLS  starts  with  initial  model  parameters  provided  by 
the  user.  These  parameters  are  then  employed  in  the  algorithm 
which  changes  these  values  by  an  iterative  method  and  evaluates 
the  resulting  fit.  If  the  fit  is  improved,  the  new  parameter  is 
accepted;  in  contrast,  if  the  fit  is  not  improved  the  previous 
parameter  is  retained  and  a  different  one  is  changed  and  the  test 
repeated  until  the  iteration  converges  to  a  criterion  of  acceptance 
evaluated  by  chi-square  (x2)  which  is  defined  according  to: 

2  _  v'  l-Zmeas  (0  —  Zmoc)e[  ( f ,  param)  | 2  ;  ^ 

x  ^  (1) 

where  Zmeas(i)  is  the  measured  impedance  at  the  frequency/],  Zm0Ciei 
{fi,  param)  is  the  calculated  impedance  based  on  the  chosen  model, 
fi  is  the  frequency  and  a,  is  the  standard  deviation  that  can  be 
equated  to  the  weight  of  the  impedance  data  points  and  this  can  be 
used  to  select  the  good  circuit  with  low  level  of  noise.  In  addition, 
the  relative  error  of  the  total  impedance  versus  frequency  indicates 
a  good  fit  and  is  highly  relevant  in  evaluating  the  acceptability  of 
the  equivalent  circuit.  The  relative  error  (normalised  to  Z)  is  defined 
as  follows: 

A|Z|(%)  =  (|Z|meas  -  |Z|calc)/|Z|meas  x  100  (2) 

where  Zmeas(i)  is  the  measured  impedance  at  the  frequency/)  and 
ZCaic  the  calculated  impedance  deriving  from  the  model  employed. 
The  overall  result  represents  an  optimisation  of  the  fit  over  the 
entire  spectrum.  Various  factors  such  as  choosing  an  incorrect 
model,  poor  estimation  of  the  initial  values  or  noise  can  prevent 
convergence  of  the  fit. 

The  data  were  also  analysed  for  Non-Linear  Least  Squares  (NLLS) 
fit  using  the  R  program.  This  was  employed  in  order  to  extract  the 
residual  standard  error  as  well  as  to  evaluate  the  fit. 

3.  Results  and  discussion 

The  Nyquist  plots  of  batteries  1  and  2  at  0%  SOC  and  15,  20,  25 
and  30  °C  are  shown  in  Figs.  2  and  3  respectively.  The  impedance 
plots  in  these  figures  may  be  represented  as  the  superposition  of 
two  impedance  arcs.  A  small  impedance  arc  (I)  can  be  seen  in  the 
medium  frequency  range  and  a  relatively  large  one  in  the  medium 
to  low  frequency  range  (II).  A  decrease  in  the  width  and  height  of 


both  impedance  arcs  can  be  observed  with  increasing  temperature, 
as  well  as  a  shift  of  the  intercepts  of  the  plots  on  the  real  axis  to 
lower  resistance  values.  The  intercept  value  decreases  from 
1.407  x  1CT3  O  at  15  °C  to  1.291  x  1CT3  O  at  30  °C  for  the  non-viable 
battery  1,  and  from  1.375  x  10~3  Q  at  15  °C  to  1.229  x  1CT3  Q  at 
30  °C  for  the  viable  battery  2. 

An  equivalent  circuit  (EC)  was  designed  for  this  study  and  is 
shown  in  Fig.  4.  A  number  of  equivalent  circuits  have  been  pro¬ 
posed  in  the  literature.  Different  combination  of  elements  gives  rise 
to  different  equivalent  circuitries  for  EIS  models  depending  on  the 
considered  chemistry  [30], 

The  fact  that  the  Nyquist  plots  in  Figs.  2  and  3  show  depressed 
semi-circles  in  the  high  to  medium  frequency  domain  suggested 
that  the  inclusion  of  CPEs  in  the  modified  Randles  circuit  would  be 
appropriate,  and  this  differentiates  our  model  from  the  DP 
approach;  our  model  may  be  considered  as  one  consisting  of 
multiple  time  constants  [31  ]  as  opposed  to  only  two.  A  process  of 
iteration  using  these  elements  at  the  anode  and  cathode  resulted  in 
the  proposed  model  assuming  a  finite  length  diffusion  component 
(Wa  and  Wc)  at  both  electrodes. 

The  proposed  EC  can  be  described  as  a  combination  of  an 
inductance  component  (L i)  connected  in  series  with  an  Ohmic 
resistance  (Ri)  which  in  turn  is  connected  in  series  with  two 
modified  Randles  circuits  comprising  a  parallel  combination  of 
constant  phase  elements  CPEa  and  CPEc  that  represent  the  double 
layer  capacitance  in  series  with  polarisation  resistances  (Rp,a  and 
RPiC),  also  known  as  charge  transfer  resistances,  and  Warburg  im¬ 
pedances  (Wa  and  Wc).  The  overall  EC  can  be  formulated  as: 
[L,  +  J?i  +  ((CPEa)/(Rp,a  +  Wa))  +  ((CPEc)/(Rp,c  +  Wc))]  where 
subscript  ‘a’  refers  to  anode  and  ‘c’  to  cathode.  Constant  phase  el¬ 
ements  (CPEs)  are  used  to  model  the  behaviour  of  imperfect  di¬ 
electrics,  and  the  electrical  impedance  of  the  CPEs  (Zcpe)  may  be 
described  by  [32]: 

ZCPE  =  l/(/w)nr  (3) 

where  w  is  the  angular  frequency  of  the  AC  signal  (rad  s_1),  and  r  is 
the  time  constant  in  F  cm~2  sn~\  with  0  <  n  <  1.  When  n  =  1,  it 
corresponds  to  a  pure  capacitor,  and  when  n  =  0  to  a  pure  resistor. 
The  Warburg  components  in  both  electrodes  represent  lithium  ion 
diffusion  processes  into  the  electrode  matrix  and  gives  rise  to  a 
complex  impedance,  Zw,  generally  defined  according  to: 

=  Rw  tanh([/wTW]n)/((/wr0)n)  (4) 

where  Rw  is  the  Warburg  resistance  (O),  rw  [33]  a  constant  in 
(seconds)  and  n  varies  between  0  and  1. 

Non-linear  Least  Squares  (NLLS)  fits  were  made  of  the  data  in 
Figs.  2  and  3,  based  on  the  EC  represented  in  Fig.  3  and  the  fitting 
parameters  are  given  in  Table  4(a)  for  battery  1  and  Table  4(b)  for 
battery  2.  The  chi-square  (x2)  values  obtained  were  all  ca.  <10~3. 
Typical  plots  of  the  error  between  the  model  and  experimental  data 
from  the  EC-Lab  software  (Zfit)  are  shown  in  Figs.  6(a— d)  and 
7(a— d).  As  can  be  seen,  the  error  is  between  0.6  and  2.4%.  More¬ 
over,  the  residual  standard  error  was  found  to  be  ca.  1CT6  using  the 
R  program,  see  Table  4(a)  and  (b). 

In  this  study,  the  high  frequency  range  with  negative  values  of 
-Im(Z)  are  ignored.  This  range  is  reported  [34]  to  be  associated 
with  an  inductive  component  (Lt)  attributable  to  the  connecting 
leads. 

The  intercept  of  the  impedance  spectrum  with  the  real  imped¬ 
ance  Re(Z)  axis  (i.e.  [-lm(Z)  =  0[)  represents  the  total  Ohmic 
resistance  [35,36]  of  the  battery  (R i)  which  is  associated  with  the 
combined  resistance  of  the  electrolyte,  current  collectors,  separa¬ 
tors,  the  electronic  resistance  of  the  active  particles  of  the  anode 
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Fig.  6.  Nyquist  plots  of  battery  1  at  0%  SOC,  at  various  temperature  and  model  simulation  ( . )  and  x2  values:  (a)  15  °C  (O)  with  x2  =  3.145  x  10-3;  (b)  20  °C  (a)  with 

X2  =  0.9034  x  10-3;  (c)  25  °C  (□)  with  x2  =  0.8595  x  10-3;  (d)  30  °C  (O)  with  x2  =  0.5227  x  10-3.  The  insert  shows  the  relative  residuals  (normalized  to  |Z|  as  a  function  of 
frequency). 
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and  cathode  materials  and  the  electrical  contacts  [37],  From  Figs.  2 
and  3,  this  Ohmic  resistance  decreases  with  increasing  tempera¬ 
ture,  which  is  most  likely  due  to  the  ionic  resistance  of  the  elec¬ 
trolyte;  which  is  temperature  dependent.  Flowever,  this  may 
include  the  bulk  resistance  of  the  electrodes. 

The  variation  of  log(l/Ri)  as  function  of  inverse  temperature  for 
batteries  1  and  2,  from  the  data  in  Figs.  1  and  2,  is  plotted  in  Fig.  5. 
As  may  be  seen  from  the  figure,  for  battery  2,  there  is  a  linear 
variation  in  agreement  with  Arrhenius  behaviour  [38],  according 
to: 

k  =  Ae-E‘!(RT'>  (5) 

where  £a  is  the  activation  energy  (J  mol-1),  R  the  gas  constant 
((J  mol-1  K-1)  and  T  the  temperature  (K).  The  slope  of  the  plot 
inFig.  5  gives  an  activation  energy  for  battery  2  of  ca.  5.3  kj  mol-1. 

This  result  is  in  reasonable  agreement  with  the  literature,  where 
such  relatively  low  values  of  activation  energy  are  associated  with 
fresh  batteries  having  undergone  few  charge/discharge  cycles.  The 
latter  generally  result  in  an  increase  in  activation  energy  due  to  an 
increase  in  the  thickness  of  the  SEI  layer  on  the  anode.  Thus  Zhang 
et  al.  [39]  have  reported  an  activation  energy  of  4.3  kj  mol-1  for  a 
fresh  battery  utilising  a  LiFeP04  cathode;  this  increased  to 
20.9  kj  mol-1  after  300  cycles.  Flowever,  an  activation  energy  of 
3.6  kj  mol-1  for  a  battery  employing  a  LiNi02  cathode  has  also  been 
reported  [40],  which  decreased  to  2.4  kj  mol-1  after  5250  cycles. 
Overall  it  appears  that  the  activation  energy  observed  is  dependent 
upon  battery  history. 

As  stated  above,  battery  1  exhibited  increased  self-discharge; 
this  coupled  with  the  non-linear  behaviour  in  Fig.  5,  suggests  that 
the  battery  had  a  fault. 

The  capacitive  arc  observed  in  the  medium  frequency  range  in 
the  Nyquisyt  plots  in  Figs.  6(a— d)  and  7(a— d),  modelled  with  a 
modified  Randles  circuit,  is  generally  attributed  to  the  SEI  layer  on 
the  anode  [41  ]  convoluted  with  diffusion  processes  (Rw).  The  sec¬ 
ond  arc,  observed  in  the  medium  frequency  region,  which  is 
generally  also  modelled  with  a  modified  Randles  circuit,  may  be 
attributed  to  the  effect  of  the  charge  transfer  resistance  at  the 
electrode  interface  [42],  also  convoluted  with  diffusion  processes 
(Rw).  The  two  capacitive  arcs  are  separated  due  to  the  different  time 
constants  for  the  processes  involved. 


The  spectra  in  Figs.  6(a— d)  and  7(a— d)  appear  as  the  super¬ 
position  of  two  impedance  arcs  in  the  medium  to  low  frequency 
range  as  discussed  above.  Both  impedance  arcs  seem  to  follow 
the  same  trend  in  that  their  widths  decrease  as  temperature 
increases.  For  battery  1,  both  impedance  arcs  merge  together  as 
the  temperature  is  increased  and  essentially  cannot  be  distin¬ 
guished  from  one  another  at  25  and  30  °C,  in  contrast  to  the 
spectra  collected  at  15  and  20  °C.  The  same  trend  is  observed  for 
battery  2;  however  the  arcs  are  only  indistinguishable  at  30  °C. 
The  separation  of  the  impedance  arcs,  and  hence  the  electro¬ 
chemical  phenomena  associated  with  them,  is  related  to  the  time 
constant  of  the  battery  [43], 

In  lithium  batteries,  the  electrochemical  processes  at  the  anode 
and  cathode  are  different  and  hence  may  be  expected  to  exhibit 
different  time  constants  which  are  obtained  from  the  following 
equation: 

r  =  (RQ o)1/n  (6) 

where  r  is  the  time  constant  in  second,  R  is  the  resistance  in  O,  Qo  is 
the  capacitance  (F). 

This  is  indeed  observed  for  both  batteries  1  and  2  as  shown  in 
Fig.  8(a)  and  (b),  where  the  logarithm  of  the  time  constant  for  the 
anode  and  cathode  of  batteries  1  and  2  are  plotted  as  a  function  of 
temperature  (at  0%  SOC). 

At  temperatures  <25  °C,  the  time  constants  associated  with  the 
processes  taking  place  in  battery  1  are  estimated  to  be  5.55  x  10-3, 
7.62  x  10-3,  6.46  x  10-4  and  8.09  x  10-6  s  at  the  cathode  at  15,  20, 
25  and  30  °C  respectively.  At  the  anode,  the  time  constants  are 
1.076  x  10-4,  0.879  x  10-4,  2.81  x  10-4  and  1.04  x  10-5  s  at  15,  20, 
25  and  30  °C,  respectively.  The  time  constants  for  battery  1  are 
smaller  at  higher  temperature  (see  Fig.  8(a)),  as  may  be  expected  for 
activation  controlled  processes.  Flowever,  it  is  clear  that  the  time 
constants  at  the  cathodes  are  greater  than  at  the  anodes  at  tem¬ 
peratures  <25  °C  in  both  batteries,  but  the  rates  become  compa¬ 
rable  at  temperatures  >25  °C. 

For  battery  2  the  time  constants  at  the  cathode  were  6.35  x  10-3, 
5.14  x  10  3,  4.34  x  10-3  and  3.76  x  10-4  s  at  15,  20,  25  and  30  °C, 
respectively  whereas  at  the  anode  the  values  obtained  were 
1.95  x  10-5,  0.967  x  10-5,  0.584  x  10-5  and  8.36  x  10-5  s  at  15,  20, 
25  and  30  °C,  respectively. 
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Fig.  7.  Nyquist  plots  of  battery  2  at  0%  SOC,  with  change  in  temperature  and  model  simulation  ( . )  and  x2  values  :(a)  15  °C  (O)  with  %2  =  1.027  x  10  3;  (b)  20  C  (A)  with 
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Fig.  7.  ( continued ). 


For  battery  2  (Fig.  8(b)),  the  time  constants  at  the  cathode  are 
also  markedly  faster  than  that  at  the  anode  at  all  temperatures 
<30  °C.  Thus  Fig.  8(a)  and  (b)  again  clearly  shows  that  the  elec¬ 
trochemical  processes  occurring  in  both  electrodes  can  be 
discriminated  for  both  batteries. 


l 

0.5  4 
0 

-0.5 

-1 

ct  -!-5 

3  -2  \ 


(b) 


-3 

-3.5 

-4 

-4.5 


Temperature  /°C 


Fig.  8.  (a)  Variation  of  log(time  constant)  with  temperature  for  battery  1  at  0%  SOC 
(temperatures  of  15,  20,  25  and  30  °C);  (b)  variation  of  log(time  constant)  with  tem¬ 
perature  for  battery  2  at  0%  SOC  (temperatures  of  15,  20,  25  and  30  °C). 


The  data  in  Fig.  8(a)  and  (b)  also  support  the  fact  that  there  is  a 
problem  with  battery  1.  Whilst  both  batteries  show  a  transition  in 
behaviour  in  terms  of  the  variation  in  time  constant  with  temper¬ 
ature,  this  transition  takes  place  at  lower  temperature  in  the  non- 
viable  battery  1. 

Battery  2  also  shows  a  transition  in  terms  of  the  time  constants, 
but  only  at  30  °C.  The  behaviour  of  both  batteries  is  reflected  in 
Fig.  9(a)  and  (b),  which  show  the  variation  in  logRpa  and  logRpc 
with  temperature.  The  resistance  of  the  anode  is  dominated  by  the 
SEI  layer  and  the  polarisation  resistance;  the  cathode  resistance  is 
totally  different  from  the  anode  [44],  Thus  the  decrease  in  the 
width  of  the  impedance  arcs  with  increase  temperature  in  Figs.  1 
and  2  is  directly  related  to  the  decrease  in  Rp  c  and  k’p  a.  In  battery 
operation  when  electrons  enter  an  electrode,  lithium  ions  diffuse 
through  the  electrode  structure  and  out  into  the  electrolyte  (charge 
transfer).  It  seems  reasonable  to  postulate  that  the  total  resistance 
of  the  anodes  in  each  battery  may  include  the  resistances  of  the  SEI 
layers  [45]  and  the  charge  transfer  resistances  [46].  However, 
experimentally  it  was  observed  that  increasing  temperature  results 
in  a  decrease  of  Rpa  for  both  batteries  (see  Tables  1  and  2)  and 
hence  to  a  decrease  in  charge  transfer  resistance  that  leads  to  the 
transition  observed  in  Fig.  9(a)  and  (b). 

The  variation  of  capacitance  with  temperature  (shown  as  log[C] 
versus  T)  is  depicted  in  Fig.  10(a)  and  (b)  for  batteries  1  and  2 
respectively.  For  battery  2,  Fig.  10(b),  the  capacitive  response  of  the 
anode  changes  little  with  temperature,  in  contrast  to  the  cathode 
which  shows  a  constant  capacitance  at  low  temperature  before 
decreasing  at  30  °C.  This  behaviour  can  be  related  to  the  charge 
transfer  taking  place  at  the  electrode/electrolyte  interfaces,  as 
capacitive  processes  are  caused  by  discharge/charge  of  the  elec¬ 
trode  surface  as  a  result  of  potential  variations.  The  change  in 
capacitance  observed  at  30  °C  may  suggest  a  structural  change  or 
SEI  formation  at  the  surface  of  the  cathode.  However  the  capacitive 
response  of  battery  1,  see  Fig.  10(a),  does  not  resemble  that 
observed  for  battery  2;  exhibiting  a  somewhat  random  variation 
with  temperature  at  both  cathode  and  anode.  It  is  clear  from 
Fig.  10(a)  and  (b)  that  we  are  able  to  discriminate  (again)  between 
the  viable  and  non-viable  batteries.  However,  further  work  at 
different  SOC,  and  an  investigation  of  self-discharge  are  required  to 
understand  and  interpret  the  data  in  Fig.  10(a)  and  (b). 

It  is  not  unreasonable  to  postulate  that  both  the  diffusion  of 
lithium  ions  through  the  electrolyte  and  the  intercalation  of  these 


702 


PS.  Attidekou  et  al.  /  Journal  of  Power  Sources  269  (2014)  694—703 


Temperature  /°C 


Temperature/°C 


3 


0.5 

0 

0  5  10  15  20  25  30  35 

Temperature/°C 


Fig.  9.  Variation  of  the  polarisation  resistance  versus  temperature  (a):  log(f?p)  at  the 
anode  (■■•  ▲  ■■■)  and  the  cathode  (  ■■  ••)  for  battery  1  at  15,  20,  25  and  30  °C  and  (b) 
log(Rp)  at  the  anode  (•■  *•••)  and  the  cathode  (  •■•  •)  for  battery  2  at  15,  20,  25  and 
30  °C. 


ions  are  facilitated  by  increasing  temperature.  Hence,  with  respect 
to  the  former  process,  we  would  expect  Rpa  and  Rpc  to  decrease 
with  increasing  temperature,  and  this  in  broad  terms,  was  indeed 
observed  in  our  data  and  by  others  [47], 

Diffusion  at  the  electrodes/electrolyte  interfaces  usually  domi¬ 
nates  the  low  frequency  domain  in  impedance  spectra,  and  this  is, 
of  course,  temperature  dependant  [48  .  The  total  impedance  of  an 
electrochemical  battery  at  a  particular  frequency  comprises  the  real 
and  imaginary  components  of  the  anode  and  the  cathode  convo¬ 
luted  together.  Persson  et  al.  [49]  have  demonstrated  extremely 
high  mobility  of  lithium  ions  between  the  layers  of  graphite.  We 
postulate  a  contribution  of  diffusion  processes  to  the  impedance  of 
the  anode  as  well  as  to  the  cathode  and  this  is  supported  by  the 
work  of  Xu  et  al.  [50]  on  a  Li/graphite  half  cell.  They  observed  a 
spike  in  the  low  frequency  domain  that  was  attributed  to  long 
range  ion  diffusion  in  the  bulk  of  the  graphite.  The  angle  with  the 
real  axis  of  the  spike  was  seen  to  change  with  cell  potential  from 
90°  progressively  to  45°,  which  is  characteristic  of  semi-infinite 
diffusion.  However,  in  order  to  determine  the  diffusional  contri¬ 
bution  of  lithium  ions  into  the  anode  and  cathode  structure, 
measured  at  low  frequency,  a  three  electrode  battery  would  be 
required,  i.e.  incorporating  a  reference  electrode  and  this  is  not  an 
option  available  for  a  sealed  Li  ion  battery.  The  only  parameter 


Fig.  10.  Variation  of  the  capacitance  versus  temperature  (a)  for  battery  1  and  (b)  for 
battery  2. 

accessible  is  the  average  diffusion  coefficient,  a  (fi  s_1),  of  the 
lithium  ions  within  the  batteries  via  the  admittance  Yo  (S  s1^2),  a 
from  which  can  be  evaluated  according  to  [51  ]: 

<t=i/(y0V2)  (7) 
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Fig.  11.  Variation  of  diffusion  coefficient  versus  temperature  at  0%  SOC,  battery  1 
(••  ■  );  and  battery  2  (•■■#•••)• 
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Fig.  10  shows  plots  of  the  a  versus  temperature  for  batteries  1 
and  2.  From  the  data  it  is  clear  that,  as  expected,  the  diffusion  co¬ 
efficient  of  both  batteries  increases  with  temperature.  In  addition, 
the  batteries  again  showed  significantly  different  response  and 
hence  could  be  distinguished  from  each  other  (Fig.  11). 

4.  Conclusions 

In  order  to  develop  a  rapid,  non-destructive,  quality  control 
method  for  lithium  battery  production,  two  commercial,  40  Ah  high 
power  lithium  ion  batteries  (viable  and  non-viable)  from  Saft  were 
investigated  using  electrochemical  impedance  spectroscopy  as  a 
function  of  temperature  at  0%  SOC.  An  equivalent  circuit-based 
model  was  formulated  to  fit  the  experimental  data  in  order  to 
extract  key  battery  parameters  for  analysis. 

The  impedance  spectra  of  both  batteries  showed  superimposed 
impedance  arcs  that  exhibited  marked  dependence  upon  temper¬ 
ature;  these  arcs  have  been  attributed  to  the  SEI  layer  and  to 
polarisation  resistance.  Because  of  the  significant  difference  in  time 
constant,  the  processes  at  the  anode  could  be  separated  from  those 
taking  place  at  the  cathode. 

In  depth  analysis  of  the  results  has  allowed  discrimination  be¬ 
tween  the  viable  and  non-viable  batteries  on  the  basis  of  a  number 
of  parameters  such  as  time  constant,  activation  energy,  polarisation 
resistances  variation  at  electrodes  with  temperature,  etc.  In  addi¬ 
tion  the  model  has  allowed  the  interpretation  of  the  data  on  the 
basis  of  physical  processes  taking  place  in  the  batteries. 

This  work  has  provided  initial  benchmark  data  on  which  a 
battery  quality  control  system  will  be  based. 
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